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ABSTRACT 

The luminosities of the centers of nearby elliptical galaxies are very low compared to 
models of thin disc accretion to their black holes at the Bondi rate, typically a few 
hundredths to a few tenths of a solar mass per year. This has motivated models of 
inefficiently-radiated accretion that invoke weak electron-ion thermal coupling, and/or 
inhibited accretion rates due to convection or outflows. Here we point out that even 
if such processes are operating, a significant fraction of the accreting gas is prevented 
from reaching the central black hole because it condenses into stars in a gravitationally 
unstable disc. Star formation occurs inside the Bondi radius (typically ~ 100 pc in 
giant ellipticals), but still relatively far from the black hole in terms of Schwarzschild 
radii. Star formation depletes and heats the gas disc, eventually leading to a marginally 
stable, but much reduced, accretion flow to the black hole. We predict the presence 
of cold (~ 100 K), dusty gas discs, containing clustered Ha emission and occasional 
type II supernovae, both resulting from the presence of massive stars. Star formation 
accounts for several features of the M87 system: a thin disc, traced by Ha emission, is 
observed on scales of about 100 pc, with features reminiscent of spiral arms and dust 
lanes; the star formation rate inferred from the intensity of Ha emission is consistent 
with the Bondi accretion rate of the system. Star formation may therefore help sup- 
press accretion onto the central engines of massive ellipticals. We also discuss some 
implications for the fueling of the Galactic center and quasars. 

Key words: accretion: accretion discs - black hole physics - galaxies: active - galax- 
ies: individual(M87) - galaxies: elliptical and lenticular - stars: formation 



1 INTRODUCTION 

Optical spectroscopy and photometry provide good evidence 
that black holes with masses of 10 s - 1O 1O M reside at the 
centers of giant elliptical galaxies (e.g. Kormendy & Rich- 
stone 1995). X-ray observations of these galaxies allow an 
inference of the interstellar gas temperature and thus an es- 
timate of the Bondi accretion radius, where the dynamics 
of the gas start to be dominated by the potential of the 
black hole. The density of the gas is also derived from these 
observations, so the Bondi accretion rate can be estimated. 
The high resolution of the Chandra X-ray Observatory has 
proved to be particularly useful for such studies. 

The bolometric luminosities of the black hole engines 
can be estimated from multi-wavelength observations, par- 
ticularly in the infrared. These luminosities are typically sev- 
eral orders of magnitude fainter than predicted by models of 



thin disc accretion at the Bondi rate in which about 10% of 
the rest mass energy is radiated (e.g. Fabian & Rees 1995). 
A particularly well-studied case is M87 (3C 274; NGC 4486) 
by Di Matteo et al. (2000, 2003) (see also Reynolds et al. 
1996a), but similar results have also been found for NGC 
1399, NGC 4472, and NGC 4636 (M49) (Di Matteo et al. 
2000; Loewenstein et al. 2001), for NGC 4697 (Di Matteo 
et al. 2000; Sarazin, Irwin & Bregman 2001) for NGC 4697, 
and for NGC 4649 (M60) (Di Matteo et al. 2000; Randall, 
Sarazin & Irwin 2004). 

In order to explain such quiescent accretors, two tem- 
perature collisionless accretion models have been invoked 
(e.g., Ichimaru 1977; Narayan & Yi 1995; Quataert & Gruzi- 
nov 1999; Narayan 2002). These models typically assume 
that the energy exchange between electrons and ions pro- 
ceeds by Coulomb collisions. The associated accretion flows 
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do not radiate efficiently if (1) the dissipation primarily 
heats the protons and (2) if the electron-ion equilibration 
time is is longer than the radial infall time. The gravita- 
tional energy is then carried by ions to the event horizon of 
the black hole, in what are known as Advection-Dominated 
Accretion Flows (ADAFs). 

For the ellipticals, minimalistic standard ADAF mod- 
els predict too high a luminosity (e.g. Di Matteo et al. 
2000). Variations of ADAFs are however still in the run- 
ning: ADAFs can be combined or replaced with a mod- 
els involving a radially dependent accretion rate and winds 
(Blandford & Begelman 1999; Quataert & Narayan 2000). 
Alternatively, ADAFs unstable to convection (called convec- 
tion dominated accretion flows - CDAFs; Narayan, Igumen- 
shchev, & Abramowicz 2000; Quataert & Gruzinov 2000) 
have the outward transport of angular momentum inhibited 
by the convection and thus their accretion rate onto the 
central engine reduced. 

There are however some unresolved issues with the 
physics of these models. Regarding CDAFs, the parame- 
ter space allowed for solutions that reduce the outward 
transport of angular momentum becomes restricted when 
magnetic fields are present and seems to require that the 
field saturates at a very low value. Otherwise the magneto- 
rotational instability dominates the outward angular mo- 
mentum transport (Balbus & Hawley 2002; Narayan et al. 
2002). It should also be noted that virtually all flavors of 
optically thin ADAFs and CDAFs require maintenance of 
two-temperature collisionless flows. The unanswered plasma 
physics questions about the viability of this situation have 
also been an important area of ongoing research (Begel- 
man & Chiueh 1988; Bitsnovatyi-Kogan & Lovelace 1997, 
Quataert 1998, Gruzinov 1998, Blackman 1999; Quataert & 
Gruzinov 2000; Blackman 2000; Pariev & Blackman 2003). 

For these reasons it is still of interest to consider other 
contributors to the quiescence of elliptical galaxy centers. 
Here we propose a very simple alternative (or companion 
process) to these accretion models. We consider the possi- 
bility that most of the accreting gas forms stars, shutting 
off the supply of fuel to the engine. In §2 we make a careful 
estimate of the Bondi accretion rate in elliptical galaxies. In 
§3 we then describe the expected properties of the disc that 
forms inside the Bondi radius, concluding that it is gravita- 
tionally unstable and should form stars. In §4 we consider 
general properties of star-forming "accretion" discs around 
black holes, describing the signatures and feedback from star 
formation, its efficiency, and the amount of gas present in 
the disc in steady state. We discuss the application of this 
model to M87, other ellipticals, the Galactic center, quasars 
and AGN in §5, and conclude in §6. 



2 THE ACCRETION RATE 

X-ray observations show that elliptical galaxies contain hot 
gas approximately at the virial temperature (e.g. Mathews 
& Brighenti 2003). Where high resolution data are avail- 
able, there is often evidence for a temperature gradient such 
that the central temperature is cooler than the outer regions 
(e.g. M87 — Di Matteo et al. 2003; NGC 4472, NGC 4636 



and NGC 5044 — Matthews & Brighenti 2003 and refer- 
ences therein). Typical temperatures on scales of several 
hundred parsecs from the galaxy centers are kT ~ 1 keV 
and so the gas is fully ionized with a mean particle mass 
of (imH with p, ~ 0.6. We define po.a = /i/0.6. This gas 
has an effective isothermal sound speed of c s = (P/p) 1 ^ 2 = 
400/i ( ^g^ 2 (fcT ( a/keV) 1 ^ 2 km s _1 , where the total pressure is 
P = nkTp and p is the density. It is possible that thermal 
pressure accounts for only a fraction, /3, of the total pres- 
sure (the remainder may include magnetic, turbulent and 
cosmic ray pressures) and we account for this by defining 
Tp = T / f3 to be an effective temperature that is a factor /3 _1 
higher than the actual gas temperature, T, i.e. as inferred 
from X-ray observations. However, to simplify the equations, 
henceforth we shall assume /3 = 1 and use Tk c v m pl ace of 
fcT/(/3keV). To allow for nonthermal pressure, one can sim- 
ply substitute Tk e v — > Tkcv//3- 

An object of mass M begins to have a significant grav- 
itational influence on the gas at a distance re, the Bondi 
radius. Note that in the case of supermassive black holes in 
giant elliptical galaxies, M may also include a contribution 
from a concentrated stellar cluster around the black hole. 
Following the notation of Shu (1992), the Bondi radius is 

_GM M 9 
r-B = — — = 26.9/i .6,oo — pc, (1) 

Cs,oc JkeV,oo 

where Mg = M/10 9 Mq and subscript oo refers to a location 
far from the black hole, i.e. at least several re- 

We can express the actual accretion rate as being some 
fraction /b of the Bondi rate 

M = /bMb = / B A 7 4 7 rp DO (GM) 2 /c 3 3i00 , 

= 0.0150/b ^n^l^MlT-^ M yr" 1 , 

where A 7 is a dimensionless factor that depends on the ef- 
fective equation of state of the gas near the Bondi radius 
(A 7 = 0.116,0.272,0.377,1.12 for 7 = 5/3,3/2,7/5,1; Shu 
1992, p437). The factor /b takes account of effects such as 
the inhibition of accretion from certain directions because 
of an outflow from the galactic central engine. 

The expected luminosity of the nuclear region inside tb 
can be expressed as 

r »> 2 o rn 1 r,43 r A 7 5/2 n^Mg _ 1 . . 

L = eMc =8.50x10 eo.i/B--^Mo.6,oc 3/2 ergs s ,(3) 

-'kcV.oo 

where eo.i = e/0.1, with this canonical value being appro- 
priate for thin disc accretion. The fact that luminosities of 
the nuclear regions of many giant ellipticals are less than the 
above value has motivated models of either inefficient radi- 
ation (small e) or inefficient accretion (small /b). We can 
relate this luminosity to the Eddington value: 

— = 6.81x10 ^fr^L-^—-. (4) 

J keV,tx> 

Now consider the heating and cooling rates at the Bondi 
radius. Their relative values determine which value of 7 is 
most appropriate in the equation of state. We assume the 
heating is only due to compression. The heating rate per 
unit volume of gas in a spherical shell at tb is 
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r = 2P B u r , B /rB = 2f^n aa kT aa P B 1 /„c Sl00 /rB 
= 1.54 x 10- 21 ./7/ u ^4T^ jOO n oo M 9 - 1 ergs cm" 3 s" 1 , [ ' 

where subscript B refers to the Bondi radius, f p = ps/Poo = 
1.44,1.55,1.64,2.45 for 7 = 5/3,3/2,7/5,1 and /„ = 
Mt,b/c s ,oo = 0.78,0.72,0.69,0.46 for 7 = 5/3,3/2,7/5,1, 
and where in the last expression we have assumed (3b = 
Poo = 1. The cooling rate for gas under conditions typical of 
the centers of ellipticals has been calculated by Sutherland 

6 Dopita (1993). For 0.0274 keV < kT < 2.74 keV and for 
metallicities that are about a factor of three greater than 
solar, the cooling is dominated by resonance line emission 
from O, Ne and Fe and can be approximated as 

A = 6.38 x 10 - 23 ^-T-°' 7 V ergs cm" 3 s" 1 . (6) 

We can use equations @ and © to evaluate when heating 
balances cooling at the Bondi radius, so that the isothermal 

7 = 1 solution is most applicable. This condition is 

w o- -, £7-2+0. 71(7-1) r ^3.21 „-3 rr,\ 

UooMq = 24.1/p' /u-MceV.oo cm , (7) 

which is also shown in Figure 1. If at a given temperature 
the combination UooMg is less than this isothermal condi- 
tion, then equations of state with 1 < 7 < 5/3 are more 
appropriate, as they account for the inability of the gas to 
cool fast enough to maintain the isothermal state. Note that 
while we have included the /3 dependence in the heating rate, 
we have not allowed for any "cooling" of the nonthermal 
component. 

Now as a specific example, consider M87. It is thought 
to harbor a black hole of mass 3.2 ± 0.9 x 10 9 M@ (this is 
the dynamical mass inside ~ 5 pc, Macchetto et al. 1997, 
but from the mass-to-light ratio stars are unlikely to con- 
tribute significantly, Harms et al. 1994), which is embed- 
ded in gas with a density near the Bondi radius of about 
0.17 ± 0.01 cm" 3 and temperature 0.80 ± 0.01 keV (Di Mat- 
teo et al. 2003). This temperature is an average of the inner 
kpc, which is a much larger scale than the Bondi radius. Us- 
ing these values the heating rate is a factor of 14.6, 4.1 times 
larger than the cooling rate for 7 = 3/2, 1, respectively (see 
also Figure 1). However the data shown by Di Matteo et al. 
(2003) show a temperature gradient so that the inner regions 
are cooler. If this trend continues all the way to the Bondi ra- 
dius then the cooling rate would balance heating, even keep- 
ing the same value of the density. Thus the existing data are 
not able to clearly distinguish whether conditions are closer 
to isothermal or adiabatic, and we must consider a range of 
possible values of 7. For the M87 system, the expected ac- 
cretion rate is O.O36/b(A 7 /O.272)M0 yr _1 and the accretion 
luminosity is 2.07 ± 1.2 x 10 44 £ .i/b(A 7 /0.272) ergss -1 . 

The observed spectrum (vF u ) from the nucleus of M87 
peaks somewhere in the vicinity of 10 pm (e.g. Di Mat- 
teo et al. 2003), but with large uncertainties. We can es- 
timate the bolometric luminosity as Lboi — iird 2 vF v ~ 
1.4 x 10 41 ergss -1 , using the observed 10.8 pm flux of 
F v = 16.7 ± 0.9 mjy (Perlman et al. 2001a) and assum- 
ing a distance of 16 Mpc (Tonry et al. 2001). Biretta, Stern, 
& Harris (1991) estimated a total core luminosity from the 
radio to the X-ray of 9.1 x 10 41 ergs s _1 . Note, however, this 
estimate is mainly based on an interpolation of flux from the 
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Figure 1. Parameter space for the balance of heating and cool- 
ing at the Bondi radius. The solid line shows the condition when 
the heating rate due to adiabatic compression and the cooling 
of 3 times solar metallicity gas are in balance, so that an accre- 
tion solution with 7 = 1 would be a reasonable approximation. 
Systems below this line are unable to cool to maintain the isother- 
mal condition and therefore 7 > 1. The reported data for M87 
are shown by the square (Di Matteo et al. 2003). A power-law 
extrapolation of the observed temperature gradient with radius 
to the location of the Bondi radius leads to the position marked 
with a circle. Given existing observational data, consideration of 
a range of equations of state, from 7 ~ 1 — 1.5 seems appropriate. 



radio to the optical bands, and is thus quite uncertain. Ap- 
plying a similar interpolation for several knots of emission 
in the jet, Biretta et al. (1991) estimate it is radiating at 
~ 2 x 10 42 ergs s _1 . Note the directly observed luminosities 
of the core and jet are only ~ 10 41 ergss -1 in each of the 
optical and X-ray bands (Biretta et al. 1991; Perlman et al. 
2001b; Marshall et al. 2002). 

The above luminosity estimates ignore kinetic energy 
that is escaping from the inner (~ 20") region. Including this 
contribution, Reynolds et al. (1996b) estimate that the total 
power of the nucleus may be about an order of magnitude 
greater, i.e. 2 x 10 43 ergss -1 , although this is sensitive to 
the 10 6 yr timescale adopted for the 5 kpc radio halo. 

Thus although the total luminosity is quite uncertain, 
both because of the limited spectral coverage for the radi- 
ated part, and the difficulty of estimating the kinetic part, it 
appears that the nucleus is underluminous by about an or- 
der of magnitude with respect to models of Bondi accretion 
with e = 0.1. 

Note that it is also possible that the nuclear activ- 
ity is varying on timescales shorter than the dynamical 
timescale at the Bondi radius, td yn = r_g/c s ,oo = 6.6 x 

iOVae.oo^s^L y r . which is ~ 3 x 1q5 X r for M87 - See 
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Waters & Zepf (2005) for a discussion on the flaring of one 
of the knots of the jet on timescales of just a few years. M87 
is the only one of the ellipticals studied by Di Matteo et al. 
(2000, 2003) that is observed to have a powerful jet at the 
present time. 

In the next section we propose a new mechanism that 
helps explain why the central black holes of massive ellipti- 
cals are underluminous: they are underfed compared to their 
Bondi rates because much of the accretion flow condenses 
into stars. Our mechanism is complementary to (not con- 
tradictory to) any additional low luminosity accretion mode 
or outflow mechanism that is operating in the immediate 
vicininty of the black hole. 



3 THE ACCRETION DISC 

The accreting material is likely to have an overall net an- 
gular momentum with respect to the central black hole. 
While the inflow is subsonic it is relatively easy for turbulent 
motions, which themselves are limited to about the sound 
speed, to transport angular momentum outwards. This has 
been observed in the numerical simulations of Abel, Bryan 
& Norman (2002) in the context of primordial star forma- 
tion from gas that is only able to cool quite slowly. Fol- 
lowing the approach of Tan & McKee (2004) we assume 
that angular momentum is conserved inside the sonic point 
r sp = vb(5 — 37)/ (47) for 7 ^ 5/3, and that the mate- 
rial passing through this point has a mean circular velocity 
^circ(r- S p) = /Kcp«Kcp(r S p). The density at the sonic point 
is p sp = poo[2/(5 - 37)] 1/(7_1) . The simulations of Abel et 
al. (2002) suggest that the circular velocities are approxi- 
mately equal to the sound speed, and adopting this as a 
fiducial value we have /k cp = (2 , y)~ 1 ^ 2 . If this is a typical 
value for the infalling gas, then the median radius at which 
material joins a rotationally supported disc is 

j-2 r sp 5-37 

r d = /Kc P r S p = — =r B ^^. (8) 

Depending on the details of the angular momentum distri- 
bution, material will be joining the disc over a range of radii. 
One possible description of this assumes uniform rotation at 
the outer boundary, which was adopted by Terebey, Shu & 
Cassen (1984) and Tan & McKee (2004). The estimate im- 
plied by eq. ||HJ is crude: the actual median radius depends 
on the details of the angular momentum distribution of the 
infalling gas. The outer radius of the disc may be several 
times larger than the median, possibly as large as the Bondi 
radius. 



3.1 Gravitational Instability and the Inevitability 
of Star Formation 

What is the fate of material that becomes centrifugally sup- 
ported at ~ Tdl For typical conditions, the densities are 
high enough so that cooling is sufficient to allow the ma- 
terial to settle into a thin disc. Local viscous processes 
and/or global instabilities can then act to transport an- 
gular momentum outwards and matter inwards. Using the 
"alpha" formalism (Shakura & Sunyaev 1973) for viscosity 



v — ac St dh = ac s d /Q, where h — c Sl d/f2 is the disc half 
thickness and Q = (GM/r 3 ) 1 ^ 2 is the orbital angular veloc- 
ity, the conservation of mass and angular momentum in a 
locally steady, thin disc imply the following relation between 
the mass accretion rate, viscosity and disc surface mass den- 
sity, E d = 2hp d , 



Md = 3nuT,d = 37ra/3dCs,d^ l S, 



(9) 



where subscript d refers to conditions in the disc and b is 
a parameter that allows viscosity to be proportional to gas 
pressure (b = 1) or total pressure (b = 0) (Goodman 2003). 
Note that use of the relation h = c Sl< j/n does not account 
for possible pressure support from radiation pressure from 
a young, massive stellar population in the disc that could 
couple to the gas via dust. The magnetorotational instability 
can yield values of a ~ 10~ 3 -10 _1 (Balbus & Hawley 1998), 
while local self-gravity can give a < 0.3 (Gammie 2001). 
If these processes are relatively inefficient then for a given 
accretion rate the disc surface density builds up to such a 
high level that the disc becomes gravitationally unstable and 
fragments into bound structures that eventually lead to the 
formation of stars. This condition can be expressed in terms 
of the dimensionless Toomre (1964) parameter 



Q 



c s ,dQ 
vrGE d 



ft 2 



2ivGpd 



(10) 



where gravitational instability occurs for Q < 1. 

Combining eqs. J5J and HOI (Goodman 2003) we have 

GM d Q = 3a/3 b d c% d (11) 

for a Keplerian disc. The numerical factor 3 becomes 2v2 
for a flat rotation curve. Although radiation pressure should 
be small compared to gas pressure at r d , justifying the as- 
sumption that f3 d — 1, stellar feedback on dusty gas could 
provide partial pressure support analogous to that of radia- 
tion pressure. Here we do not consider this further and set 
Pd = 1. 

First consider a disc that is heated purely by viscous 
dissipation and is cooled by radiation from an optically thick 
surface. The effective temperature is 



3 

8ITCT 



4.63 



3 

kcV, 



( A-r V ( &~, 2 
I 0.272 J \^5-37 



1/8 



• \ 1/4 
J 

3/4 

K, 



(12) 



where a is the Stefan-Boltzmann constant and we have used 
Eq. 101. If the disc is optically thick (an assumption that we 
check later) then the midplane temperature is given by 



T 4 



(13) 



assuming the heating rate per unit mass and the opacity per 
unit mass are constant. As seen from eq. 11211 we are likely 
to be in a regime where dust grains can form and survive 
(i.e. Td < 2000 K). The Rosseland mean opacity per unit gas 
mass then depends on T d as well as the gas-to-dust ratio and 
the properties of the grains. For example at 100 K and for 
a gas-to-dust ratio similar to the Milky Way's the opacity is 
about 1 cm 2 g" 1 (Li & Draine 2001). 

The value of the disc surface mass density is 
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153a, 



-4/5 



xT, 



B "°° 0.272 
3/5 



3/5 



9/10 M 4/5 



-3/20 -1/5 4/5 / 8-v 2 V 7 ( J_\ 



Mo. 6 

-3/5 



(14) 



g cm 



where fj,d is the mean particle mass in the disc in units of 
mH. The midplane temperature is 

2/5 



T d = 12.7 



a 0.3 



0.272 



1/10 ^27/80 



m"'° & 



9/10 (15) 

K. 



Note that this is the minimum temperature that an opti- 
cally thick disc can have, since the calculation assumes only 
viscous heating. Including the factors that depend on 7, the 
fiducial temperatures are 321, 78 K for 7 = 3/2, 1. Such a 
disc is indeed optically thick (jd = S^Kd > 1), if dust can 
form that has properties similar to typical Milky Way dust 
grains. 

What is the value of Q for the optically thick disc? 
Expressing equation (I I U in terms of the fiducial parameters 
of Bondi accretion in a typical giant elliptical galaxy, we have 

-2/5 



: 4.86 x 10-V d /5 ^ 3/4 °- 7/1 ° 



v ,," 47 /20 M -17/10 T 291/160 
X Mo.6,oc IV1 9 1 keV,oo 



r-z/5 ( X y \ 
0.3 JB \^0.272 I 



-2/5 

"'OC 

27/20 



(16) 



The disc is likely to be unstable with respect to its self grav- 
ity: for example applying the parameters appropriate to M87 
yields Q = (0.600, 0.0175)alQ f~ 2/5 K,f /40 for 7 = 3/2,1. 
Note the relatively large value of Q in the 7 = 3/2 case is 
a consequence of the small value of r^. If a radius closer to 
the observed disc size in M87 (~ 100 pc) is adopted then 
Q< 1. 

In summary, the accretion of gas to supermassive black 
holes in the centers of elliptical galaxies is likely to lead 
to the formation of a gravitationally unstable gas disc. The 
natural expectation is that this disc will fragment into bound 
objects: stars. We consider this process and its consequences 
in the next section. 



4 STAR-FORMING DISCS AROUND 
GALACTIC BLACK HOLES 

We expect that gravitational instability in the disc leads 
to fragmentation and star formation. This process removes 
mass from the gaseous disc, reducing £<£. The newly formed 
stars heat, and thus stabilize, the remaining gas. 

First consider the reduction in the surface density of gas 
due to star formation. If the stars remain in a thin disc, the 
system will still be prone to gravitational instabilities. These 
will heat the stellar population relative to the gas since the 
latter can cool. As the scaleheight of the stellar distribution 
becomes larger than that of the gas, the stars have less and 
less influence on the gravitational stability of the gas. Fur- 
thermore, for continued star formation, the gas itself must 
be self-gravitating, not just the star-gas system. This effect 
can be seen in the star-forming discs of typical disc galaxies: 
the threshold for star formation is still adequately described 
by the Toomre parameter, Q, that depends on the total sur- 
face density of gas (Martin & Kennicutt 2001). The presence 



of the more massive stellar disc raises the particular value 
of Q below which star formation occurs, but only by fac- 
tors of a few or less (Jog & Solomon 1984). One way that 
this can occur is via the local concentration of gas in stellar 
spiral arms (e.g. Kuno et al. 1995). In summary, as T, d is 
reduced via star formation, the gas disc becomes more and 
more stable with respect to its self-gravity. 

An additional effect of the reduction of £d is a smaller 
optical depth for the disc's cooling radiation. Cooling be- 
comes most effective when Td ~ 1. Other things being equal, 
we would then expect a somewhat cooler disc midplane tem- 
perature than in the optically thick limit, given above. How- 
ever, now the disc is being heated by stellar photospheres 
that are much hotter than Td- For a uniform distribution of 
gas over the face of the disc, the optical depth through the 
disc to this radiation from dust absorption is still large, even 
if the optical depth to thermal radiation at T d is ~ 1. 

The amount of heating per mass of new stars depends 
on their initial mass function (IMF) . The upper part of the 
IMF appears to approximately follow a Salpeter (1955) dis- 
tribution, djV/dm oc m" 2 ' 35 , for 1 M© < m < 100M Q with 
a dearth of stars beyond the upper mass limit. The IMFs in 
the Galaxy (Salpeter 1955; Kroupa 2002) and in metal-rich 
starbursts (Schaerer et al. 2000) appear to be similar in this 
regime. Below a solar mass the Galactic IMF flattens some- 
what (e.g. Muench et al. 2002), but much of the total mass 
is still contained in sub-solar-mass stars. There are no ob- 
servational constraints on this part of the IMF in starburst 
galaxies. 

For the above IMF, while most of the mass is in low- 
mass stars, most of the energy injection into the ISM comes 
from high-mass stars. However the value of the high-mass 
cutoff in the IMF does not significantly affect the overall 
feedback properties of the population, so long as it is greater 
than ~ 6OM0 or so. The uncertainty in the low-mass end in- 
troduces an overall uncertainty in the normalization of feed- 
back per unit stellar mass of factors of a few: for example if 
the Salpeter IMF extends down to 0.1 Af© then the feedback 
per stellar mass is reduced by a factor of 0.39 from the case 
with a lower-mass cutoff of 1 Mq . 

Goodman & Tan (2004) have argued that gravitational 
instability in the inner regions of quasar accretion discs (at 
radii of a fraction of a parsec) may lead to very massive stars. 
Conditions there are very different from those at a typical 
value of rd of the Bondi-fed discs considered in the present 
paper. In particular, the inner quasar discs are closer to be- 
ing dominated by radiation pressure as opposed to gas pres- 
sure; temperatures are too high to allow the presence of dust 
grains; and these regions are only marginally self-gravitating 
(by definition). We expect that these differences are impor- 
tant enough to lead to very different IMFs. In short, given 
existing data we feel that adopting a Salpeter IMF is the 
most reasonable approximation to make for star formation 
in Bondi-fed gas discs of ellipticals. 

For a Salpeter IMF from 0.1 to 100 Mq, the lu- 
minosity associated with a given star formation rate, 
M», is L* ~ 3.9 x 1O 9 (M*/M yr" 1 )!,© = 1.5 x 
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1O 43 (M»/M yr _1 ) ergss -1 (Leitherer et al. 1999). 1 After 
about 30 Myr, this stellar population produces a super- 
nova at a rate of 7.8 x 10~ 3 '(M* /Mq yr _1 )yr _1 , causing 
a mean injection rate of mechanical energy of ~ 2.3 x 
10 41 (M,/M yr 1 ) ergs s 1 . The mechanical luminosity of 
stellar winds is estimated to be about a factor of 6 smaller. 
Thus the maximum energy input from stars is about e* = 
E,/(M,c 2 ) = 2.6 x 10" 4 . The number of H-ionizing photons 
produced is (9.2, 7.9, 7.2) x 1Q 52 (M»/M yr~ 1 )s~ 1 for metal- 
licity Z = 1.0,2.0,3.0^0 (Smith, Norris, & Crowther 2002), 
where the highest metallicity value is based on a power law 
extrapolation from the lower two values. Assuming case B re- 
combination in the ionized regions with no escape of ionizing 
photons yields an Ha luminosity of (Kennicutt, Tamblyn, & 
Congdon 1994) 



0.99 x 10 41 /tra P ,i(M./M Q yr" 1 ) ergs s" 



(17) 



with the normalization being for the Z — 3Zq case and 
where /trap.i is the fraction of ionizing photons that are 
trapped in the disc's HII region and do not escape into the 
hot interstellar medium of the galaxy. This estimate ignores 
dust absorption in the disc, which would lower the numerical 
factor in the above equation. Note that a recombination flow 
of one solar mass per year would produce an Ha luminosity 
that is a factor of 4.4 x 10 -4 times smaller than the amount 
produced by one solar mass per year of star formation (with 
the above Salpeter IMF). 

The Ha flux from the circumnuclear disc of M87 is 
4.4 ± 1.5 x 10" 14 ergs cm -2 s" 1 (Ford et al. 1994) (see Figure 
2), corresponding to a luminosity of 1.35 x 10 39 ergs s _1 for 
an assumed distance of 16 Mpc (Tonry et al. 2001). This im- 
plies a star formation rate of 1.4 x 10 / t ~ ; Mq yr -1 using 
the normalization of equation 11711 . Note that this estimate 
is a lower limit because it does not allow for any internal 
extinction in the star-forming disc that would attenuate the 
observed Ha flux. 

If we assume that the bolometric luminosity of the cir- 
cumnuclear disc in M87 is similar to the observed peak in 
i/Lv, i.e. ~ 1.4 x 10 41 ergss -1 (Perlman et al. 2001a; see 
§2), and is due entirely to star formation, then the above 
relations imply a star formation rate of 9.5 x 10 -3 Mq yr -1 . 
Figure 3, described in more detail in the next section, shows 
how the bolometric luminosity due to star formation com- 
pares to the observed spectrum of the M87 nucleus. 

Despite the uncertainties in the IMF, internal extinc- 
tion and equation of state of the accreting gas, the key re- 
sult from the above star formation rate estimates is that 
they are both roughly consistent with the inferred Bondi ac- 
cretion rate. This agreement suggests that the star-forming 
accretion flow model provides a plausible mechanism for pre- 
venting significant gas accretion onto the central black hole, 
at least in M87. This is in addition to the fact that a rela- 
tively thin disc is observed to be present on scales just inside 
the Bondi radius (see further discussion in §5.1). 

1 This is actually the luminosity of a stellar population that has 
been continuously forming stars at l.OM yr _1 for the last 10 7 yr, 
which is approximately the age at which ionizing feedback reaches 
a constant level. 



- 




i 



1" 



Figure 2. Adapted from Ford et al. (1994): HST F658N on- 
band image minus HST F547M off-band image of the M87 nucleus 
designed to highlight Ha emission. North is towards the top and 
east to the left; the famous synchrotron jet (not visible here) 
extends to the NW. The black hole is at the center of the image 
and is surrounded by a disc of material: the velocity difference 
from 1" to the NE to 1" to the SW is 1000 km s" 1 (Harms 
et al. 1994; Macchetto et al. 1997). The presence of this disc, 
including its spiral structure and dust lanes, as well as its total 
Ha luminosity, is evidence in support of a model of circumnuclear 
star formation. 



4.1 Estimate of the Star Formation Efficiency and 
Residual Accretion to the Black Hole 

Here we present a simple one zone model that describes the 
star formation efficiency of a Q = 1 disc as a function of 
the coupling of star formation feedback to the gas. From 
equation l|lljl we have 



Q = 



3a 
GM d 



kT d 
lid 



3/2 



(18) 



For an optically thick disc, T d ~ 3/t<jST^/8 and its lu- 
minosity is L ~ 2ltr d aTe S = /trap.L-k. ^ /trap,L3.9 X 

1O 9 (M*/M yr _1 )L , assuming a negligible contribution 
from viscous heating, valid at r ~ r d if we are to achieve 
a quasi-stable disc with Q ~ 1. In the last step above, we 
have used the IMF-dependent relation between luminosity 
and star formation rate discussed previously. Now by using 
the above relations, including eq. (1 lOt . assuming that Q = 1 
and that most of the Bondi accretion rate goes into stars, i.e. 
M* ~ M, we derive the mass accretion rate in the disc that 
can be supported against self-gravity by stellar feedback: 
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M d = 2.27 x 



2.686 x KTW: 



(iv^) 3/7 (Tofe)" 3/2 ^0^-; 

/ K ti/trap,L-^l 



M 9M o.6,oo 0.272 
3/2 



3/7 (19) 



where Li = L»/3.9 x lO 9 !/© is the luminosity generated by 
a star formation rate of 1 Mq yr -1 for our adopted IMF. 
This is the residual amount of gas accretion that can be 
stabilized with respect to its own self-gravity by transform- 
ing most of the initial mass flux into stars. We have written 
eq. 119H in two stages so that different constraints can be 
used: either the observed disc size and the Bondi accretion 
rate to the disc (or any general mass feeding rate to the 
disc), or the properties of the hot gas in the elliptical (in 
this case the theoretical estimate for the disc size is used; 
eq.0. For M87, in the first case the accretion rate in the disc 
is (7.0 x f0" 7 ,1.3 x lO- 6 )ao.3(Kd/trap,L/s/irf 4 ) 3/7 M yr" 1 
for 7 = 3/2, f with rd = 100 pc, while in the second case it 

is (1.36 xlO" 4 , 9.25 xlO- 6 )a . 3 («:d/tr a p,L/BMd 4 ) 3/7M oy r " 1 
for 7 = 3/2, 1. These are much lower than the initial Bondi 

accretion rates of (0.036, 0.148)/ B Mq yr -1 . 

If the disc is optically thin to its cooling radiation, then 
- ^eff,d/ r d = TeB,d/{ K d^d), and the equation for the 
mass accretion rate in the disc becomes 



M d = 8.57 x 10- 5 a . 3 

1/3 



r.? N l/6 



(^)"(T^)~ 1/6 M yr- 

/ f r t 2 w \ 1/3 \ ZU ) 

2.63 x 10- 5 ao. 3 f W^^ D .,. B .«, ^ 



Now for M87, in the first formula of J2lJ the 
accretion rate in the disc is (2.3 x 10 _5 ,3.7 x 
10- 5 )a .3K7 tr ap,L/ S Md 4 ) 1/3 M© yr" 1 for 7 = 3/2,1 
with Td — 100 pc, while in the second formula it is 
(4.2 x HT 5 ,4.7x 10- 5 )ao.3(Kd/trap,L./BMd 4 ) 3/7 M Q yr" 1 for 
7 = 3/2,1. 

We can judge which of the optically thin or thick lim- 
its is appropriate by the value of Ed = Mdfl/(37raCg d ) = 
(Md/(3a)) 1/3 fi/(7r(GQ) 2/3 ) and the disc temperature, 
which sets Kd- For the optically thick case the surface density 
is 



7.05 x 10" 



v — k — 



Mq yr 



1/7 



V 100 pc) 

while in the optically thin case it is 



g cm 



0.0236 



A/, 



1/9 



pr) gem . 



(21) 



(22) 



For a dust-to-gas ratio three times that of the Milky Way's, 
Av — 600E so that the above two numerical factors corre- 
spond to Av = 4.2, 14, respectively. The effective tempera- 
ture at the disc surface is 



r off ,d - {^^) - 25.8 ^ M0yr -i ) (23) 

x(lofe) 1/2 K 

For most elliptical galaxies, the disc will be optically thin 
to its cooling radiation in the outer parts, and will make a 
transition to the optically thick regime in the inner region. 
Note that the discs are expected to be quite optically thick 
to their stellar radiation field, although clustering of young 
stars and ISM inhomogeneities will act to reduce this. Be- 
cause of these complications we have left the factor /trap.L 
as a free parameter, but expect it to have a value of order 
unity. 

Despite the uncertainties in the parameters and the de- 
pendence on the adopted disc size, the estimates suggest 
that only a small fraction of the initial mass accretion flux 
through the disc would be supported against star forma- 
tion and accreted by the central black hole. The amount of 
the reduction implied by equations l|19|l and (1201 is more 
than sufficient to explain the low luminosity of M87's nu- 
cleus. However, the residual accretion rates predicted by the 
above simple model are so small that other processes that 
have not been considered, such as mass supply from stel- 
lar winds, may now be relatively important. Another caveat 
is that the above estimates only consider thermal pressure 
support in the gas disc and ignore the additional stabilizing 
effect of magnetic fields. 

Note that some material will be expelled from stars 
as they undergo and complete their stellar evolution. For 
Salpeter-like IMFs, most of this matter will be returned af- 
ter a relatively long timescale: for our adopted Salpeter IMF 
from 0.1 to 100 Mq, only 14% of the mass is in stars with 
masses greater than 8 Mq. Irrespective of the timescale of 
mass return, most of this material would be recycled back 
into the gas disc rather than accreted directly by the black 
hole. This is because most stars are forming and evolving at 
distances from the black hole that are still quite large com- 
pared to the sizes of their stellar wind bubbles or planetary 
nebulae, and the specific angular momentum of the returned 
gas is the same as that of the stars. 

Other authors have highlighted the issue of self-gravity 
in the accretion discs that feed supermassive black holes, 
mainly in the context of quasars that are accreting at much 
higher rates, approaching their Eddington limits (Schlosman 
& Begelman 1989; Goodman 2003; Sirko & Goodman 2003; 
Thompson et al. 2005). Goodman (2003) pointed out that 
stellar heating of the disc was unlikely to be able to stabilize 
the disc, because the energy release per mass of stars formed 
is relatively low. We have extended this line of reasoning to 
conclude that most of the mass flux goes into stars, which 
are then able to stabilize a small residual fraction of the gas 
in an accretion disc to the black hole. The implications of 
this model for the fueling of quasars are discussed in §5.4. 

4.2 Estimate of the total mass of gas in the disc 

One prediction of our model of (Bondi-fed) star-forming ac- 
cretion flows is the presence of cold, presumably molecular 
gas. To estimate the amount, we assume that a steady state 
has been achieved so that the star formation rate in the disc 
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balances the Bondi accretion rate to the disc. The disc mass 
is Md — Tvr^'Sd- In the optically thick case this is 

1/7 



M d = 1.06 x 10 b 



Kd/tlap.L-klMg 



M, 



l4 



M , (24) 



(25) 



M e yr" 1 

while in the optically thin case it is 
M d - 3.55x10 ^ ^ m yr-i 

\ 100 pc) 

For M87 the disc mass is then (1.3, 1.6) x 
lO^Kd/trap.LilMd 4 ) 177 M e (7 = 3/2,1) in the opti- 
cally thick case and (4.1, 4.8) x 10 8 (K^ 1 /t ra p,L£i/v') 1/9 M Q 
(7 = 3/2, 1) in the optically thin case. 

As an alternative method, we use an empirical rela- 
tion between star formation rate and disc gas mass to esti- 
mate the properties of the gas disc. For a sample of about a 
hundred galactic and circumnuclear discs, Kennicutt (1998) 
found an empirical linear relation between the surface den- 
sity of star formation and the product of the total mean gas 
surface density and the mean orbital angular velocity: 

Ssfr ~ 0.017E d fi, (26) 

where Esfr = M*/(7rr^). It is not clear if this relation ap- 
plies to Bondi-fed circumnuclear discs in ellipticals as the 
orbital timescales are much shorter than in the sample of 
Kennicutt and the rotation curves may be closer to Keple- 
rian (e.g. M87 - Macchetto et al. 1997), rather than being 
flat or rising. This may lead to somewhat different star for- 
mation efficiencies, at least in the context of shear- dependent 
models of star formation (Tan 2000). 

Nevertheless, using eq. 1261 we can estimate the surface 
density of the gas disc: 



M e yr- 1 



M, 



-1/2 



100 pc 



-1/2 



g cm 



-2 



(27) 



In the case of M87 this takes values (3.7, 15.3) x 
10~ 3 gcnr 2 ee (18,73) M e pc" 2 for 7 = 3/2,1, which is 
in the middle of the range of values covered by Kennicutt's 
sample. The total gas mass in this disc is 



2.77 x 10' 



Af, 



M d = 
1/2 



, 0.017H / oa \ 
r d \3/2 „, r (28) 

M yr- 1 9 V 100 pc 

which takes values (0.6,2.3 x 10 6 ) M Q in the case of M87 
with 7 = 3/2, 1. The coincidence of the two different meth- 
ods of estimating the disc mass, suggests that this feedback 
model may have relevance to setting the normalization of 
the empirical Kennicutt (1998) relation. 

Note that if the disc were disrupted by a particular 
event, such as an outburst from the AGN, then it would 
take ~ 10 7 — 10 8 yr to reform, if supplied at the Bondi rate. 



5 DISCUSSION 

We have presented a model of star-forming accretion flows 
with particular attention to the flow inside the Bondi radius 
of the central black hole of an elliptical galaxy. Our approach 



was as follows: first we calculated the accretion rate, showing 
the explicit dependence on the effective equation of state of 
the gas. We then estimated the heating and cooling rates of 
the gas in the accretion flow near the Bondi radius, showing 
that idealized accretion solutions 7 ~ 1 — 1.5 are probably 
quite a reasonable description of the actual accretion flow, 
at least in the case of M87. We then argued that angular 
momentum conservation inside the sonic point of the flow 
leads to formation of a disc on scales that are a factor of 
a few or so smaller than the Bondi radius. The disc was 
shown to be gravitationally unstable if the heating is due to 
viscous processes in the disc. We found that self-gravity is 
so strong that most of the mass flux will be converted into 
stars. The energy input from these star can stabilize only a 
small residual fraction of the initial mass flux: i.e. the accre- 
tion rate to the black hole is substantially reduced below the 
Bondi accretion rate. We argue that this is a major cause 
of the "low luminosity" of black holes observed in giant el- 
liptical galaxies, where "low luminosity" is in reference to 
a model of thin disc accretion at the Bondi rate as deter- 
mined from the conditions in the kpc-scale, X-ray-emitting 
gas. A radiatively-inefficient accretion flow may exist at the 
center of the disc, much closer to the black hole. However, 
its inefficiency must be argued in the context of an accretion 
rate and a luminosity determined at a much smaller scale, 
for example where the disc is no longer self-gravitating. 

The predictions of this model include the presence of a 
star-forming gas disc (i.e. containing relatively cold molecu- 
lar gas) around the black holes of ellipticals and the presence 
of hot, young stars in these regions. These may reveal them- 
selves by spectral features (e.g. clustered Ha emission from 
gas ionized by massive stars; stellar atmosphere/ wind fea- 
tures, e.g. from Wolf-Rayet stars; or stellar ultraviolet con- 
tinuum, although this may be swamped by AGN emission) 
or by the occurrence of type II supernovae. For example a 
typical star formation rate at the center of a giant ellipti- 
cal galaxy is expected to be ~ 0.1 M© yr -1 , about a factor 
of 30 less than the star formation rate of a Milky- Way-like 
disc galaxy. For the standard Salpeter IMF adopted here, 
the type II supernova rate from such an elliptical galaxy is 
~8x 10" 4 yr" 1 . 

Nayakshin (2004) considered how the presence of many 
young stars, left over from a previous quasar phase, can 
affect the accretion and appearance of low-luminosity AGN. 
However, he did not consider the possibility that the stars 
are currently forming in these systems. 

We now summarize the potential application of the star- 
forming accretion disc model to M87, other giant ellipticals, 
the Galactic Center, quasars and other AGN. 



5.1 Application to M87 

We have used M87 as an example to illustrate the model 
throughout the text. Here we summarize this application, 
discuss the existing evidence that supports the model, and 
suggest future observations that can confirm or refute its 
validity. 

The central black hole of M87, with mass (3.2 ± 0.9) x 
10 9 Me (Macchetto et al. 1997), should have a Bondi accre- 
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tion rate of 0.036 - 0.15 M© yr" 1 (this range corresponds to 
uncertainties in 7 = 3/2 — 1) from the hot (0.8 keV), low- 
density (n e — 0.17 cm -3 ) gas (Di Matteo et al. 2003) at the 
center of the giant elliptical (uncertainties in the black hole 
mass and gas temperature could also change the accretion 
rate by factors of three or so). The Bondi radius is about 
110 pc, but could be twice as large if the high end of the 
black hole mass determination is valid and the temperature 
in this inner region is actually 0.5 keV (see discussion in §2). 
The expected median disc radius is several times smaller 
than the Bondi radius: our fiducial case has it smaller by 
factors of 36 and 4 for 7 = 3/2, 1. The actual outer disc 
radius would be larger than the median by factors of a few, 
and in any case the exact size depends on the details of the 
initial angular momentum distribution of the hot accreting 
gas, which is not well constrained. Since the Ha disc seen by 
HST is ~ lOOpc in extent (Ford et al. 1994; see Fig. 2), this 
is probably the best estimate of the disc scale in this system. 
Note also that the HST observations show some filamentary 
structure out to several hundred pc, so that the accretion 
process may not be quite a simple as given by the purely 
Bondi solution. 

Despite the caveats, the coincidence of the observed disc 
size with the estimated Bondi radius supports a model in 
which that disc is forming from a Bondi accretion flow with 
initial circular motions of order the hot gas sound speed. 

The Ha disc was described by Ford et al. (1994) as hav- 
ing "spiral structure" and "dust lanes" (see Fig. 2), which 
supports the theoretical conclusion that it is gravitationally 
unstable (§3) and contains cool gas that is liable to undergo 
star formation (§§3,4). The total Ha luminosity implies a 
star formation rate (assuming a standard Salpeter IMF) 
within factors of a few of the Bondi accretion rate at which 
the disc is being fed. This is consistent with the implied star 
formation rate being equal to the Bondi accretion rate given 
the uncertainties. 

We can model the bolometric luminosity associated 
with this star formation rate. Figure 3 shows the observed in- 
frared and radio spectrum of the M87 nucleus, much of which 
results from nonthermal emission from the AGN (see §2). 
Care must be taken since the fluxes at different wavelengths 
often correspond to regions of different scales. Also shown in 
Figure 3 is a simple black body spectrum of the gas and dust 
disc heated by star formation at rates 0.036 — 0.15 A/© yr~\ 
corresponding to luminosities (1.6 — 5.9) x 10 s L©. The tem- 
perature is set by assuming this luminosity is emitted over 
the two faces of a disc of radius 100 pc. This model is highly 
idealized, since it assumes a single temperature, ignores ad- 
ditional heating from the interstellar radiation field of the 
galaxy, ignores deviations from a pure black body spectrum, 
such as PAH emission features in the range 1 — 20 /Ltm, 
and assumes all the light from stars is reprocessed by dust. 
The overall luminosity normalization is also IMF-dependent. 
Nevertheless it is interesting that these simple models are 
now in a position to be tested by high resolution infrared and 
sub-mm observations that are now feasible with the Spitzer 
Space Telescope and the Sub-millimeter Array (SMA). 

The presence of cold, presumably molecular gas in the 
disc with total mass ~ 10 6 Mq could be searched for with 
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Figure 3. Continuum spectrum of the nucleus of M87 shown as 
F v (top) and vF v (bottom). Data in order of increasing wave- 
length: HST optical (Ho 1999); Gemini 10 /Urn (Perlman et al. 
2001a); Circles — SMA calibration data (sensitive to flux on 
scales of the order 2"); Triangles — IRAM 30m data (sensitive 
to flux on scales of order 11" and 28" at 1.3 mm and 3 mm, re- 
spectively (data also extracted from the SMA calibrator list); 
Crosses — VLBI data at 100 GHz (Baath et al. 1992) and 22 GHz 
(Spencer & Junor 1986) sensitive to flux on scales of 0.2" and 
~0. 00015". The solid curve is a black body spectrum from a 
star-forming disc with the star formation rate equal to the up- 
per limit of the Bondi feeding rate of the disc (0.15M© yr — 1 ) and 
a standard Salpeter mass function down to 0.1 M©. The tem- 
perature used here is 16 K set by having the emitting area be 
equal to that of the observed Ho disc with radius ~ 100 pc. 
This is a lower limit as it neglects additional heating from the 
evolved stellar population of the galaxy. The dotted curve is the 
equivalent model for the lower limit of the Bondi accretion rate 
(0.036 Mgyr" 1 ). 

mm interferometers. Such a search is particularly well suited 
to the capabilities of the SMA and ALMA. 

To estimate the efficiency of the central engine we need 
to know the black hole's accretion rate. As noted in §2, to 
explain the observations requires a reduction in the accretion 
rate by a factor of ~ 10 from the Bondi value (this factor 
is itself uncertain by about an order of magnitude), or a 
similar reduction in the efficiency of energy liberation from 
the fiducial value of e — 0.1. Note that most of the energy 
is thought to be released as non-radiating kinetic energy 
of the jet/outflow. The reduction in black hole accretion 
rate from the Bondi rate due to star formation is also quite 
uncertain. We estimated (eq. I20H a possible reduction by 
factors of up to ~ 1000. However, as discussed in §4, this 
simple estimate is based on a simple one-zone model for 
the disc and does not include effects such as accretion of 
some very low angular momentum gas from the larger-scale 
galaxy and the return of material evolved stars. In summary, 
the observations of the M87 nucleus can be explained by 
models of efficient (e ~ 0.1) accretion if a large fraction 
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(~ 90%) of the initial mass flux to the disc forms stars. 
While such a high fraction is certainly allowed by the simple 
estimates of self-regulated star formation we have presented 
in §4, quantitative estimates of the residual accretion rate 
are difficult. 



5.2 Other Elliptical Galaxies 

Optical imaging with HST has been used to detect dusty 
discs in elliptical galaxies. Jaffe et al. (1993) reported a disc 
of cool gas and dust with radial extent of 65 pc around the 
nucleus of NGC 4261 (3C 270), another elliptical galaxy in 
the Virgo cluster (see also Martel et al. 2000 and van Bem- 
mel et al. 2004). Bower et al. (1997) reported a disc, with 
radius 82 pc, of ionized gas in M84 (NGC 4374 = 3C 272.1). 
Jaffe et al. (1994) surveyed 14 Virgo ellipticals, finding clear 
evidence for dust discs (seen in absorption) in several. Note 
this method of detection could be sensitive to the orien- 
tation of the disc and other properties of the nucleus: for 
example in this study no dust disc was seen in absorption in 
the case of M87. They note that essentially all their systems 
exhibit a change in morphology on scales of ~ 100 — 200 pc 
from the nucleus: either in the extent of an absorbing dust 
disc or a change in the brightness profile. We identify this 
scale with the Bondi radius, inside which a cool, star-forming 
disc can form. Van Dokkum & Franx (1995) found that dust 
"lanes" were present in the majority (~ 80%) of a sample 
of 64 early-type galaxies. These were usually perpendicular 
to the axis of any radio jet present in the system. Carollo 
et al. (1997) and Tomita et al. (2000) also found evidence 
for nuclear dust absorption in the majority of the galaxies 
they analyzed. Martel et al. (2004) found nuclear discs of 
dust and ionized gas in a large fraction of the nine early 
type galaxies they studied with HST- ACS. Again the typ- 
ical scales of these features are ~ 100 pc. They estimated 
star formation rates of (0.6 — 6) x 10 -3 A'/oyr -1 assuming all 
the Ha emission is produced by star formation and using the 
normalization of Kennicutt (1998). Lauer et al. (2005) found 
dust in the centers of about half their sample of 77 galaxies. 
Only one of their "dusty" galaxies had dust that was more 
prominent on scales greater than 4"from the center than 
closer to the center; most of the dust is highly concentrated, 
and in all the systems where the dust is clearly in a disc 
or a ring, no dust clouds are visible external to the outer 
edge of the disc or ring. Significantly, no dust was seen in 
relatively low luminosity, low mass ellipticals (My > —21), 
which might be expected since the Bondi accretion rate in- 
creases with galaxy mass, assuming the black hole mass is a 
fixed fraction of the galaxy mass. 

Jaffe & McNamara (1994) detected radio absorption in 
HI and CO from neutral gas in the disc of NGC 4261 and 
estimated a mass ~2x 10 5 Mg, within an order of magni- 
tude of the theoretical mass estimates we have made for the 
disc in M87. Knapp & Rupen (1996) surveyed 42 elliptical 
galaxies for CO emission, concluding that the detection rate 
was just under 50%. However, since the half-power beam 
of these observations was about 30", these detections most 
likely do not relate to molecular gas forming inside the Bondi 
radius, but to material formed from evolved stars or asso- 



ciated with the infall of gas-rich satellite galaxies. Young 
(2002) detected CO emission in 5 out of 7 ellipticals, but 
again, because of the resolution of the observations, the fo- 
cus was on ~kpc scales. Vila-Vilaro, Cepa, & Butner (2003) 
detected CO(3-2) emission from 6 out of 10 early-type galax- 
ies that already showed CO(2-l) emission, with a beam of 
about 20". They inferred typical star formation rates (from 
Ha) of a few tenths of a solar mass per year, but again this 
is mostly on scales larger than the Bondi radius (nuclear Ha 
emission was not considered in this study). 

In summary, there is much evidence that elliptical 
galaxies contain dusty gas discs on scales approximately 
equal to their Bondi radii, i.e. ~ 100 pc. It is less clear if 
star formation is actively occurring in these discs. One could 
make a more quantitative test of whether central dust con- 
centrations are related to the Bondi radius by comparing the 
size of the dust features with the vb- If the gas temperature 
near the Bondi radius is proportional to the virial tempera- 
ture, then Too oc T v \ r oc al oc Afg a p, where <r* is the galaxy's 
stellar velocity dispersion and M ga i is its mass, and the last 
step uses the empirical relation of Faber et al. (1997). As- 
suming a constant black hole to galaxy mass fraction, i.e. 
M oc Mgai , the Bondi radius then scales as re oc Mg a f 8 . 

Long term star formation at the scale of the Bondi ra- 
dius (~ 100 pc) will have implications for the morphologies 
of the cores of elliptical galaxies. For example it may account 
for the ~ 100 pc scale discs seen in some early type galaxies 
(e.g. van den Bosch, Jaffe, & van der Marel 1998; Krajnovic 
& Jaffe 2004). It may also account for the flattening in the 
radial profiles of "core" elliptical galaxies (e.g., Lauer et al. 
2005): this occurs on scales of ~ 100 pc in galaxies that tend 
to be the most massive ellipticals. 



5.3 The Galactic Center 

The black hole at the Galactic center has a mass of ~ 
3 — 4 x 10 6 Mq, appears to be embedded in hot gas (as 
seen in X-rays), yet is underluminous by several orders of 
magnitude with respect to thin disc accretion at the Bondi 
rate (Baganoff et al. 2003, and references therein). Some, 
but probably not the majority, of the diffuse X-ray emission 
on the arcsecond scales relevant to Bondi accretion may be 
the dust-scattered halo of the X-ray emission from inside 
the black hole's Bondi radius (Tan & Draine 2004). Basic 
ADAF models have been ruled out for the Galactic center 
due to polarization measurements which indicate the den- 
sity is even lower than that predicted by ADAFs (Bower et 
al. 2003). 

This system is probably more complicated than the 
supermassive black holes in elliptical galaxies because the 
Bondi radius is only ~ 0.1 pc in size, which is relatively 
small compared to the feedback scales associated with stel- 
lar wind bubbles and supernova remnants from individual 
massive stars. For example, the dominant source of mass in 
the region is likely to be from stellar winds from young, mas- 
sive stars. Some of these are even present inside the Bondi 
radius (as estimated for the diffuse X-ray emitting gas) and 
because of the relatively low mass of the black hole, these 
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winds with speeds ~1000 km/s may be unbound and act to 
prevent diffuse gas from accreting (Quataert 2004). 

A further complication is that there are large amounts 
of neutral gas in the vicinity (Herrnstein & Ho 2004; Christo- 
pher et al. 2005). The sporadic infall of this material may 
well dominate the time-averaged mass accretion rate. 



5.4 Quasars and AGN 

The issue of self-gravity and quasar accretion discs has 
been discussed by Schlosman & Begelman (1989), Good- 
man (2003) and Thompson et al. (2005). Our model herein 
has been developed with low-luminosity galactic nuclei in 
mind, for which the central regions are being fed at the 
Bondi rate from their host elliptical galaxies. The accre- 
tion rates to the black holes powering quasars approach 
the Eddington limited value: Md.E = AuGMIe/ {ne.s.ce) ~ 
22(Ze/£o.i)A/9 M0yr _1 , where n e . s ~ 0.4cm 2 g" 1 is the elec- 
tron scattering opacity, Ie = L/Le is the ratio of the lumi- 
nosity to the Eddington value, and e = L/(Md<?) = O.leo.i 
is the radiative efficiency of the accretion disc. 

To have these accretion rates be the residual mass flux 
from a star-formation supported disc on 100 pc scales re- 
quires star formation rates that are ~ 10 6 times larger 
(eq. IHJ, which are unrealistic. However, this estimate as- 
sumes the star formation time is short compared to the ra- 
dial advection time in the disc. Thompson et al. (2005) show 
that the opposite may be true in quasars, in which case the 
initial outer-disc mass feeding rate (or equivalently the star 
formation rate) need only be factors of ~100 greater. 

There may be additional complications, such as the sup- 
ply of some very low angular momentum material to the 
inner disc, which is not gravitationally unstable. 

The problem of the mass supply to quasars requires 
further work. The possibility remains that self-gravity in a 
disc is so severe that black holes need to be supplied at much 
smaller radii, perhaps through stellar disruption events (e.g. 
Hills 1975; Goodman 2003). 

In lower luminosity AGN, such as Seyferts and LIN- 
ERs, there is much evidence for nuclear star formation (e.g. 
Terlevich & Melnick 1985). These accretion discs are also 
expected to be self-gravitating, and there is some observa- 
tional evidence to support this (e.g. Kondratko, Greenhill, 
& Moran 2005). While the fueling of these nuclei may have 
little to do with Bondi accretion (because of the presence of 
large amounts of cold gas in an extended disc), the model of 
an accretion disc supported by stellar feedback (§4) probably 
is relevant, predicting a general correlation of star formation 
and AGN activity. 



6 CONCLUSIONS 

Observations of the central regions of elliptical galaxies sug- 
gest that supermassive black holes that reside there are fed 
at the Bondi rate, of a few hundredths to a few tenths of a 
solar mass per year. For standard thin disc accretion, these 
black holes are very underluminous given this mass accretion 



rate. Previously, alternative accretion models invoking inef- 
ficient radiation of gravitational energy, radially-dependent 
accretion rates, and outflows have been proposed as solu- 
tions to this problem. In this paper we have suggested an 
additional process, likely occurring, which can be of com- 
parable or potentially greater importance: Bondi-fed star- 
forming discs. Accretion occurs at the Bondi rate, but be- 
cause the disc that inevitably forms inside the Bondi radius 
is very gravitationally unstable, most of the mass flux reach- 
ing it turns into stars. Stellar feedback is relatively weak so 
that only a small residual amount of gas can be stabilized to 
reach the black hole at the center of the disc. Our estimates 
of this residual mass flux are uncertain, but appear to be 
potentially small enough that standard thin disc accretion 
models can avoid the constraints imposed by the observed 
low luminosities. ADAFs, CDAFs, or other models with out- 
flows may still apply in the very inner regions, but with a 
smaller outer boundary value for their inward mass flux. 

Our model is highly simplified: we mostly consider only 
a single typical scale for the disc radius, whereas in real- 
ity gas will reach the disc at a range of radii depending 
on the initial angular momentum distribution; our quan- 
titative estimates depend on the assumed stellar IMF; we 
have adopted the standard alpha approximation for viscous 
stresses in the disc; and we have mostly considered only ther- 
mal pressure support in the disc, ignoring magnetic pressure. 

As a result of these simplifications, it is vital to com- 
pare the model as closely as possible to real systems, the 
best observed being M87. We have described how the ob- 
servational data on the nucleus of this galaxy are consistent 
with the possibility that more than half of the Bondi accre- 
tion flow is depleted into stars: a thin disc of gas and dust, 
apparently self-gravitating, is seen inside the Bondi radius 
and the inferred star formation rate from the Ha emission 
is consistent with the Bondi accretion rate. This conclusion 
is not contradicted even when current estimates of the me- 
chanical luminosity of the M87 jet are incorporated into the 
residual central engine accretion rates. Additional observa- 
tions to detect molecular gas scales inside the Bondi radius 
are now feasible with the SMA. This instrument is also ca- 
pable of detecting continuum emission from dust heated by 
the expected star formation. 
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